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A discrete charge transfer in a small tunnel junction where Coulomb interactions are important 
can excite electron-hole pairs near the Fermi level. We use a simple model to study the associated 
, nonequilibrium properties and found two novel effects: (i) for junctions with electrodes of the same 

■ electronic properties, a leakage current exists within the Coulomb gap even when the environmental 

Q\ \ impedance is infinite; (ii) for junctions with electrodes of different electronic properties, the differ- 

ential conductance diverges when a net interaction between conduction electrons is attractive, and 
it is strongly suppressed for a net repulsive interaction. 

(J^ ' PACS numbers: 75.fO.Jm, 75.fO.Lp, 75.30.Ds. 
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Small tunnel junctions show unique features associated 
with a discrete charge transfer across the energy barrier 
0. These systems are being studied in detail, and have 
striking implications for other fields of physics [||. 
f***- ■ In this Letter, we focus on tunneling processes occur- 
ring at energies below the elementary charging energy 
Ec — e 2 /2C, e and C being the electronic charge and 
the junction capacitance. This is the energy range within 
which small tunnel junctions show their most interesting 
properties, such as Coulomb blockade. Our main results 
in this regime are: 

i) Junctions whose electrodes have the same electronic 
properties, show a leakage current within the Coulomb 

O ■ gap even when the environmental impedance is infinite. 
This leakage current smoothes the current-voltage (I-V) 
characteristics and reduces the Coulomb offset extrapo- 
lated back from the high-voltage part of the I-V curve. 

ii) In junctions with electrodes of different electronic 
properties, the I-V characteristics show rectifying fea- 
tures within the Coulomb gap, and under certain circum- 
stances, conductivity diverges at zero voltage, exhibiting 
a superconducting-like behavior. 

While most discussions are restricted to single junc- 
tions, this Letter also presents extensions to multijunc- 
tion arrays. 

A distinguishing feature of our analysis is that it takes 
account of the fact that electron tunneling in small junc- 
tions strongly perturbs the electronic states close to the 
Fermi level. The conventional theory assumes that the 
junction is in a new equilibrium state immediately after 
a tunneling event. This implies that the electronic wave- 
functions which describe the electrons near the Fermi 
level are not affected by the change in the charge state 
of the junction, and only a rigid shift in the energy lev- 
els must be considered. The validity of this assumption, 
in fact, depends on the inhomogeneities in the potential 
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generated after the charge build-up, over scales compa- 
rable to the mean free path of the electrons. When the 
inhomogeneities cannot be neglected, the dependence of 
the electronic wavefunctions on the charge state of the 
junction must be included in the calculations. The im- 
portance of this effect increases as the electrodes become 
smaller. 

If the change in the electronic wavefunctions is not neg- 
ligible, each charge transfer is followed by sizable shake- 
up processes ||. The Nozieres-Dominicis problem de- 
scribes, in a similar way, the ejection of a core electron 
which changes the charge state of an atom within a metal 
A substantial shakeup of the conduction electrons 
follows, leading to singularities in the photoemission and 
X-ray absorption spectra. As in the junction, the sys- 
tem experiences a fast change -the modification of the 
charge state- and a resulting slow relaxation. The main 
alteration in the junction case is that the transferred elec- 
tron that 'triggers' a time-dependent perturbation is it- 
self part of the Fermi liquid. This difference, however, 
does not qualitatively change the above picture because 
the final state interaction is caused, not by the tunneling 
electron itself, but the change in the collective coordinate 
which describes the charge state. 

For more precision, we separate the tunneling event 
into two processes having very different time scales: (i) 
an electron near the Fermi level in one electrode tun- 
nels to the other, giving rise to a charge build-up, and 
(ii) low-energy electron-hole pairs are excited in response 
to this quick change in the potential. This idea can be 
implemented mathematically by separating the degrees 
of freedom of the junction into two classes: a collective 
one, which lumps together all 'fast' variables, and the 
dynamics of the electrons at the Fermi level, which can 
be described in terms of independent electron- hole pairs. 
The dynamics of the fast variable are associated with 
changes in the charge state of the junction. We define 
this collective variable as Q and its conjugate variable 
as 4> such that [<j>, Q] = ie, where e is the electronic 
charge. The wavefunction of the junction can be writ- 
ten as a superposition of states specified by the variable 
4> and the occupancy of the levels around the Fermi en- 
ergy: \4>, {nk} >■ Charge relaxation occurs over energy 
scales comparable to the plasmon energy, ~ leV, while 
we use the quasiparticle description for processes at en- 
ergies, or voltages, comparable to the charging energy of 
the system, Ec ~ 1 — lOmeV, so that our distinction is 
justified. 

Including all these effects, the junction can be de- 
scribed by the Hamiltonian ||, 

Ti-c + Hl + H-r + T~Ct + 

(Q - Qxf 
2C 

^2<^k,R,Lct,R,L c k,R,L 
k 



n = 

H c = 
T~Lr,l = 
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H T = ie i ' l> J2 c i,R c k',L + h.c. 

k,k> 

T-C V = Vr(Q - Q x ) ^ C k,R C k'-,R 
k,k' 

-Vl(Q - Qx) / ] c k.L C k',L (!) 
k,k' 

Here Q x stands for any offset charge which may exist 
at equilibrium. The differences between Hamiltonian (1) 
and more conventional treatments are the separation of 
slow and fast degrees of freedom, and the introduction 
of the shakeup processes through the last term || . Note 
that the tunneling hamiltonian TLt is defined such that 
modification in the charge state is accompanied by a 
change in the distribution of quasiparticles at the Fermi 
level. 

We can eliminate the slow degrees of freedom, leaving 
only the charge. The influence of the electrons on the 
evolution of the charge can be expressed in terms of elec- 
tronic propagators . The inclusion of the shake-up pro- 
cesses implies that the propagators have to be computed 
for the nonequilibrium ground state of charged electrodes. 
We will use only their asymptotic behavior at long times 

§: 

G R ' L (t) = -i^ll(E c t) g ^ (2) 

In equation (2), the effects of the potential on the elec- 
trons at the Fermi level, due to vr and vl, are described 
by the phase shifts induced at the Fermi level, <5r ; l, and 
g R X = 25r L /it - {Srl/tt) 2 . 

The contribution that has the anomalous exponent 
arises from the shake-up. An effective coupling between 
tunneling processes separated by time t is described by 
the product G R (t) x G L (t), which depends on a single 
parameter, g — —gR — gh- This approach, which is based 
on the definition of an appropriate tunneling density of 
states that includes the effects of the interaction, is rather 
general, and can be applied in many contexts 0. 

We now consider the case of a single junction with 
electrodes of the same electronic properties. After a tun- 
neling event, the potentials acting on the states at the 
Fermi level are equal in magnitude but opposite in sign, 
and so are the phase shifts, i.e., Sr — —5l = 5. Hence, 
we find that g — 2(S/ir) 2 > 0. In the small coupling 
limit, t/Ec <C 1, the current operator, ed(NR — Ni)/dt, 
is 

j = iete l * 4,R c k' ,l + h - c - (3) 

k,k' 

Using Kubo's formula, the conductivity at finite frequen- 
cies and zero temperature is given by 

= s i < Q|jy > i 2 (4) 
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where the summation is taken over all excited states of 
energy u>. 

The numerator in (4) is the Fourier transform of a 
correlation function of the type < j(t)j(t r ) >. We now 
express this correlation function as a series expansion in 
powers of Ht- Each insertion of this operator generates 
a kink in the time evolution of the charge. A schematic 
representation of a given term in this series is shown in 
Fig. 1. 

We can use the standard analogy with a one dimen- 
sional system with logarithmic interactions Q . The con- 
tribution of the path shown in Fig. 1 then reduces to the 
interaction effects between two charge dipoles at separa- 
tion t. It decays over a long time like t(~ 4 ~s\ Taking 
the Fourier transform, and inserting it in (4), we obtain 

„M ~ it*N R (E F fN L (E F f (Jt) 2+ * 

where Nr(Ef) and Nl(Ef) denote the respective den- 
sity of states at the Fermi level in the right and left 
electrodes. This formula is valid at low frequencies, 
huj < E c . 

The existence of a finite conductance at frequencies 
within the Coulomb gap also implies a finite differential 
resistance in the same region. Equation (5) is a measure 
of the probability that the junction absorbs energy in 
quanta given by tku. Thus, at low voltages, we expect 
I cx V 3+ ^, as shown in Fig. 2. 

Our results imply that a single normal tunnel junc- 
tion is always dissipative. The decay of the 'Coulomb 
blockade state' cannot be significantly reduced by low- 
ering the temperature or increasing the resistance in the 
external circuit. We stress that, in an experiment de- 
signed to measure the FV characteristics, this current 
within the Coulomb gap does not violate the conserva- 
tion law of energy. To determine the I-V characteristics, 
the junction has to be attached to a reservoir in which 
dissipation can take place. Otherwise, even in the ab- 
sence of Coulomb blockade, electrons at the electrodes 
cannot be described by Fermi-Dirac distributions. When 
the circuit is attached to a macroscopic battery, processes 
which do not change the electrochemical potential differ- 
ence V across the junction are not forbidden because, 
as in the Mossbauer effect, the battery instantaneously 
compensates the work eV. This situation is similar to 
what we have encountered in the low-impedance envi- 
ronment. In this case, spontaneous fluctuations in the 
junction charge due to the electromagnetic environment 
induces tunneling within the Coulomb gap; this process, 
however, does not violate the conservation law of energy 
because the battery again compensates the energy neces- 
sary for tunneling In the high-impedance case, how- 
ever, the second-order process in tunneling is inhibited, 
and therefore our fourth-order process comes to play a 
major role. 



4 



The same result can be expressed in terms of the 
fluctuation-dissipation theorem. To make the chemical 
potentials at the two electrodes unequal, the system has 
to be brought out of equilibrium. The perturbation re- 
quired to do so has an intrinsic time scale, given by Ti/eV 
(a voltage can be induced by an ac magnetic field, for 
instance). Processes which have longer relaxation times 
play no role in the steady state reached by the system. 
That the absence of these processes give rise to a voltage- 
or temperature-dependent renormalization of the junc- 
tion capacitance has been discussed extensively. By the 
same analysis, these processes can also contribute to the 
imaginary part of the response function, leading to the 
dissipation we have shown above. 

The existence of a finite subgap conductance will also 
reduce the Coulomb offset extrapolated back from the 
high-voltage part of the I-V curve. Most of the current 
will come from tunneling processes which involve energies 
between Eq and the applied voltage. They will give rise 
to a linear dependence of I on V. A current due to the 
small subgap conductance also contributes, reducing the 
Coulomb offset. Both the smoothing of the I-V charac- 
teristic and the reduction in the Coulomb offset seem to 
be in agreement with experimental observations [jio] |l2| . 
Finally, dissipation will be accompanied by heating ef- 
fects ||. 

The analysis presented here can be readily extended 
to many junction arrays. The frequency (or voltage) de- 
pendence of the conductivity is not changed because it 
is simply related to the interaction between two charge 
dipoles in a manner similar to that depicted in Fig. 1. 
The prefactor, on the other hand, has to be replaced by 
(h/e 2 R) N . This term increases exponentially with the 
number of junctions for high resistance junction arrays. 
It will compete with other processes already discussed 
in the literature, like macroscopic quantum tunneling of 
charge |l4| - |l6f| . Physically, the latter implies two coher- 
ent tunneling events at two different junctions, and is 
drastically different from the shakeup we discuss here. 
It seems an interesting question to separate these two 
sources of subgap conductance. From the requirement 
of coherence, we conclude that macroscopic tunneling of 
charge should be very sensitive to changes in the con- 
ductance of a single junction. Shakeup processes, on the 
other hand, should always be present, independently of 
the ratio between the different conductances. 

We next consider junctions with electrodes of different 
electronic properties for which 5r Sl- We concentrate 
on the high conductance limit, a — h/e 2 R ^> 1. In this 
case, it is convenient to start with the effective action: 

b e ff = ,L± a ± T c 

Jo dT Jo dT l T - T ^M ail ^\ + 2 ™% ( 6 ) 

where the parameter r c is the cutoff which is initially set 
equal to the charging energy. Since the two possible elec- 
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tron jumps between the electrodes are described by dif- 
ferent propagators, they are associated with the two sign 
choices in the action. The first term is complex, indicat- 
ing an instability of the 'neutral' junction. A real par- 
tition function is recovered by assuming a temperature- 
dependent offset charge Q x at the electrodes. We have 
defined two different conductances, a±. At high volt- 
ages or temperatures, we have a+ = a_ = ao, where 
a = h/(e 2 R). 

The high frequency modes in the phase can be inte- 
grated out in (6), leading to the scaling equations |l7| 

gg± = | 2 
dln(eV) a± 7T 2 

a-(g R + gL) (7) 



dln(eV) 



Here, the second equation is needed to keep the parti- 
tion function real at all scales. The first equation is ob- 
tained by scaling separately the contributions from the 
two types of hopping processes. 

These equations have a nontrivial fixed point at a cr it = 
2/{-K 2 \g R + g L \). For initial conductances above this 
value, the conductances within the Coulomb gap scale 
towards infinity as the scale of frequencies (or voltages) 
is reduced. Below this threshold, we recover the scaling 
equations for symmetric tunnel junctions already known 
||. Integrating (7), we obtain 

/ V \±(.gn.+9L) 

a±(V) = a ( f^-J «o > «crit 

a±(V) = a - ^ log (fe-) a -Ca crit (8) 

The first of these equations leads to asymmetric I-V 
characteristics. The differential resistance at zero van- 
ishes for the positive sign and diverges for the negative 
sign. The qualitative features of the I-V characteristics 
in the two cases are shown in Fig. 3. 

The scaling equations (7) are valid in the limit of high 
conductances. This condition is always satisfied when 
the frequency (or voltage) is scaled towards lower values 
in the high-conductance side of the I-V characteristics. 
This effective 'normal-state superconductivity' is a novel 
effect which merits further theoretical and experimental 
study. 

The range of parameters required for its observation 
can be inferred from the previous analysis. The differ- 
ence in phase shifts between the two electrodes must be 
large. This implies that they must be very different, and 
that the potentials induced by their respective charge 
distributions must be local in space. A negative phase 
shift indicates that an additional charge pulls states be- 
low the Fermi level, while a positive shift has the opposite 
effect. One of the electrodes should therefore show a net 
repulsion between charges, while the other should have 
a net attraction. Similarly, it can be shown that a di- 
vergent conductivity arises from tunneling between two 
attractive Luttinger liquids H] . Finally, the normal state 
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conductance, ao, has to be large, meaning good contact 
between the electrodes, or many channels through which 
electrons can tunnel. The latter possibility, however, will 
tend to reduce the phase shifts by spreading the potential 
of a unit charge over the different channels. We therefore 
think that this effect would be better observed in small 
tunnel junctions with good connections. 
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FIG. 1. Path integral formulation of the calculation of the 
< j(t')j{t) > correlation function (see text). 
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FIG. 2. I-V characteristics of a normal, symmetric tunnel 
junction. 

FIG. 3. I-V characteristics of an asymmetric tunnel junc- 
tion, where V c = Ec/eexp(— 7r 2 (a — 1)/2) is the renormalized 
offset voltage. Solid line: ao <ti a C rit- Broken line: ao 2?> a cr it. 
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